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Neutron scattering experiments have been performed on Sri4Cu24 04i which consists of both 
chains and ladders of copper ions. We observed that the magnetic excitations from the Cu02 
chain have two branches and that both branches are weakly dispersive along the a and c axes. 
The ui-Q dispersion relation as well as the intensities can be reasonably described by a random 
phase approximation with intradimer coupling between next-nearest-neighbor copper spins J=ll 
meV, interdimer coupling along the c axis J c =0.75 meV, and interdimer coupling along the a axis 
J a =0.75 meV. The dimer configuration indicates a quasi-two-dimensional hole ordering, resulting 
in an ordering of magnetic Cu 2+ with spin-i and nonmagnetic Cu, which forms the Zhang- Rice 
singlet. We have also studied the effect of Ca substitution for Sr on the dimer and the hole ordering. 
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I. INTRODUCTION 

Sri4Cu2404i consists of both two-leg ladders of copper ions and simple CuO2rCbfins0i as shown in Fig. 1. Numerous 
experiments showed that the two-leg ladder has an excitation gap of ~35 meV,tTTj which is expected theoretically^ An 
important feature of this compound is that stoichiometric Sri4Cu2404i contains hole carriers. It has been reported 
that most of the holes are localized in the chain and some exist in the ladder .ErlHl When Sr 2+ sites are substituted by 
Ca + ions, total number of the holes in the sample is unchanged but holc&iii the chain are transferred from the chain 
to the ladderE£Eil and the system shows an insulator-to-metal tran^ioroi 2 ]. Superconductivity was also observed in 
Sr .4Cai3. 6 Cu 2 4O4i below T C =10K under a high pressure of 3 GPalI3. 

In this paper we are only concerned with the magnetic properties in the chains. Figure 1(a) shows the Cu02 chains 
in this compound. The copper ions are coupled by almost 90° Cu-O-Cu bond along the c axis. Each chain is well 
isolated from each other along the a and b axes. As-jpcntioned above, there are localized holes in the chain. It-is 
expected that the hole spins are localized at oxygenE2l and couple with copper spins to form Zhang-Rice singleto. 
The nonmagnetic Cu sites play an important role to form a dimerized state in the chain. First, the dimer is formed 
between next-nearest-neighb©j-|-Cu 2+ spins along the c axis. The exchange interaction (~10 meV) is mediated via 
a nonmagnetiCj-ZR singlet .BtLalia Then the question is how the dimers are arranged and interact with each other. 
Matsuda et al.tB interpreted that each dimer is separated by one ZR singlet along the c axis (Model I) as shown in 
Fig. 1(c). They also interpreted that the excitations at low-Q {Lchain <0.4) originate from a dimerized state in the 
chain and those at high-Q {Lchain >0.6) from a dimerized state due to the interladder coupling since two branches 
were observed for the gap excitations at low-Q but only one branch at high-XjL. 

NMR studies revealed the microscopic properties of the dimerized stateJj'EJ It was reported that both magnetic 
Cu 2+ and nonmagnetic ZR singlet exist in the chain. It was also found that a NMR peak originating from ZR singlet 
gradually splits into two peaks below ~200 K. The results were discussed with two dimer -models. One is Model I 
mentioned above. This model was supported by Cox et al. by their x-ray measurements Jl3 In another model the 
dimers, which are formed between next-nearest-neighbor Cu 2+ spins, are separated by two nonmagnetic ZR singlets 
along the c axis (ModeJ II) as shown in Fig. 1(d). This model is consistent with the results of the electron diffraction 
study by Hiroi et al&B in which the Cu0 2 chains show a modulated structure with five times larger unit cell along 
the c axis in stoichiometric Sri4Cu2*Q«. However, neither models cannot completely explain the NMR results. ■— ■ 

Recently, two unpublished paperal3'E3 on neutron scattering experiments came to our attention. Ecclcston et al$3 
explained the modulation of the excitations ancL-uitcnsities along the c axis convincingly by a simple model of an 
alternating chain with weak interdimer couplingEU {J%j J\ ~ —0.1, where J\ and J2 are intradimer and interdimer 
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coupling, respectively), which is equivalent to Model II. It is noted that two parallel branches along the c axis are not 
separated in the experiments because of insufficient resolution. Regnault et alEB briefly reported the measurement of 
two parallel branches in a wide range of Q (0.05< L c h a in <0.925) along the c axis, which suggested a weak, coupling 
along the a direction. The specific coupling along the a direction was recently reported by Cox et alS3 in their 
synchrotron x-ray study of the charge ordering at low temperatures. 

Since dimerized states sometimes occur in spin (S) h one-dimensional Heisenberg antiferromagnets such as the 
spin-Peierls state, it is quite natural to assume that the dimerized state in the chain of Sri4Cu2404i is also caused by 
a quantum effect in S=^ one-dimensional Heisenberg antiferromagnet. On the other hand, almost isolated magnetic 
dimers were found in CaCuGe206EjEj even though the geometrical arrangement of the magnetic moments seems 
three-dimensional. It was also claimed that magnetic moments wiuch are geometrically closest does not necessarily 
couple most dominantly in VODPO4 • ^D20ca and (VO)2P207l3 because of the strong superexchange pathways 
through a covalently bonded PO4 group. These results suggest that unexpected pathways could become important 
in order to completely understand the magnetic properties of magnetic materials. I— . 

We have searched for a simple model of dimers with weak couplings along both c and-4 directions.Ea Somewhat 
surprisingly, a specific combination of the two couplings, as described by Leuenberger et aZ.EZIfor CssCi^Brg, produces 
simple and elegant neutron scattering cross sections which describe properly the measured dispersion and intensities. 
In this paper, we have studied uj-Q dispersion relation perpendicular to the chain direction in considerable detail. 
By applying a random phase approximation (RPA) treatment, we found that the interdimer coupling along the a 
axis (J a =0.75 meV) is also importanto as well as the interdimer coupling along the c axis (</ c =0.75 meV). Model 
cross section will be given after the experimental data are presented. The dimer arrangement along the c axis is 
well described with Model II. The dimer configuration we found in this study indicates a quasi-two-dimensional hole 
ordering, resulting in an ordering of Cu 2+ and ZR singletin the ac plane. Ca substitution for Sr makes the magnetic 
excitation peaks broader as in the case of Y substitutions probably because the hole ordering is very sensitive to the 
hole number and the long-range dimer formation becomes disturbed. 

II. EXPERIMENTAL DETAILS 

The single crystals of Sri4_ ;E Ca 2 ;Cu2404i (x—0 and 3) were grown using a traveling solvent floating zone (TSFZ) 
method at 3 bars oxygen atmosphere. The dimension of the cylindrically shaped crystals is about 5 x 5 x 30 mm 3 . 
The effective mosaic of the single crystal is less than 0.4° with the spectrometer condition as described below. The 
Sri4Cu2404i crystal is the same one as used in Ref. It is expected that Sr and Ca are distributed homogeneously 
in SrnCa3Cu2404i since the lattice constants systematically change and the linewidth of the nuclear Bragg peaks 
does not change when the ratio of Sr and Ca is changed. The lattice constants of Sri4Cu2404i and SrnCa3Cu2404i 
are et=11.472 A and c=27.551 A and a=II.430 A and c=27.487 A at 15 K, respectively. The lattice constants are 
consistent with those obtained with powder samplesB 

The neutron scattering experiments were carried out on the ISSP-PONTA spectrometer installed at the 5G beam 
hole of the Japan Research Reactor 3M (JRR-3M) at Japan Atomic Energy Research Institute (JAERI). The horizontal 
collimator sequences were 40'-40'-S-80'-80'. The final neutron energy was fixed at Ef=14.7 meV. Pyrolytic graphite 
(002) was used as monochromator and analyzer. Contamination from higher-order beam was effectively eliminated 
using pyrolytic graphite filters after the sample. The single crystals were mounted in a closed cycle refregirator and 
were oriented in the (h, 0, 1) scattering plane. As described in Ref. [I], there are three different values for the lattice 
constant c (c un i VC rsai=10 x c c h a in=7 x ci a dder)- Since we will show the magnetic and structural properties in the chain, 
Cchain will be used to express Miller indices. 

III. EXPERIMENTAL RESULTS 
A. Magnetic excitations from Sri4Cu2404i chain 

1. Low temperature dispersion 

Figure 2 shows the typical neutron inelastic spectra at (H,0,L) in Sri4Cu2404i measured at 15 K. One or two 
distinct excitation peaks am^pbserved in the (H, 0, L) scattering plane as in the (0, K, L) scattering plane which was 
previously reported in Re£aM The solid lines at (2,0,-0.20), (2.75,0,-0.20), (2,0,-0.75), and (3,0,-0.75) are 
fits to two Gaussians and those at (2, 0, —0.70) and (2.75, 0, —0.70) are the fits to a single Gaussian. The excitation 
peak positions are changed when H or L is changed, meaning a dispersion along the a and c axes. The width of the 
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excitation peaks varies at different Q positions purely due to the resolution focussing effect. At (3,0, —0.75), where 
resolution focussing is almost perfect, the peak width has a resolution limited value of ~1 meV full width at half 
maximum. The observed excitation energies and energy-integrated intensities at (2,0, L) and (3,0,L) (0.1< L <0.85) 
are plotted in Figs. 3 and 4, respectively. The closed circles represent the data from the magnetic peaks with lower 
excitation energy and the open circles from the ones with higher excitation energy. The solid and broken lines are 
the results of model calculations, which will be presented in the next section. Both at (2,0, L) and (3,0,-L) the 
dispersion relations-are similar and has a_periodicity of 0.2 r.l.u. along the c axis which is consistent with the results 
by Eccleston et alL2 and Regnault et alcn Actually the dispersion relation is almost identical to the one observed in 
the (0, K, L) zoneJlj The most interesting feature is that the intensities and the ratio of the two peaks are changed 
with H as well as L. The intensities of the excitation peaks with lower excitation energy are more intense than those 
with higher excitation energy at (2,0, L) (L <0.25). On the other hand, the intensities of the excitation peaks with 
higher excitation energy are more intense than those with lower excitation energy at (2, 0, L) (L >0.6). The intensities 
at (2, 0, L) are almost identical to those observed in the (0, K, L) zone.t3 The intensities at (3, 0, L) show completely 
opposite behavior from those at (2, 0, L). It is noted that the intensities in Ref. [if] is the sum of intensities from the 
two excitations. 

Figures 5 and 6 show the observed excitation energies and energy-integrated intensities at (H, 0, —0.2) and 
(H, 0,-0.7) (2< H <4), respectively. The closed circles represent the data from the magnetic peaks with lower 
excitation energy and the open circles from the ones with higher excitation energy. The solid and broken lines are 
the results of model calculations, which will be presented in the next section. There also exists a dispersion along the 
a axis with a periodicity of 2 r.l.u. The two dispersion curves along the a axis do not run parallel as ones along the c 
axis but cross at H = | and J. The intensities of the excitation peaks with lower excitation energy are larger than 
those with higher excitation energy in a range of 2< H <2.5 and 3.5< H <4 at (H, 0, —0.2). On the other hand, in a 
range of 2.5< H <3.5 the intensities of the excitation peaks with higher excitation energy are larger than those with 
lower excitation energy. The intensities at (H, 0, —0.7) show completely opposite behavior from those at (H, 0, —0.2). 
Thus, both the dispersion and the intensities depend on H and L, suggesting that there are nonnegligible magnetic 
correlations along the a axis in addition to those along the c axis. 



2. Model Hamiltonian 



In order to analyze the observed dispersion relation and energy-integrated intensities we used a model Hamiltonian 
for the dimers formed between next-nearest-neighbor Cu 2+ spins which are weakly coupled along the a and c axes as 
shown in Fig. 7. The Hamiltonian involves three Hcisenbcrg Hamiltonians with intradimer coupling J, interdimer 
coupling in the same chain J c , and interdimer coupling between the adjacent chains J a : 

<i,j> <i,j>' <i,j>" 

Here, < i,j > and < i,j >' are nearest-neighbor and next-nearest-neighbor spins in the same chain, respectively. 
< i, j >" is nearest-neighbor spins between the adjacent chains. This type of diieters, namely weakly, coupled dimers, 
can be well described with RPA. In fact, the magnetic excitations in, CssC^BrgKZl and BaCuSi20gE3 are successfully 
described with RPA treatment. The dispersion relation is given bj£3 

^acoustic/optic^) = [J2 + J . R ^ (JJ^) ± J a \ T<j ( g ) | )] 1/2 ; 

mT , i - ex p(- J/T) 

fi(T) = n -n 1 = i + 3exp( _ J/T) , 

r c (g) = 2cos(2tt/), 

T a (q) = 2cos(tt/i) ■ exp(i27rW), (2) 

where no, n\, and S are thermal populations of the singlet ground state and of the first excited triplet and the distance 
along the a axis between nearest-neighbor spins at the adjacent chains as shown in Fig. 7. The solid and broken 
lines in Figs. 3(a)-6(a) represent the calculated dispersion relation of the acoustic and optic modes, respectively. 
Equation (2) with J=ll mcV, J a =0.75 meV, and J c =0.75 meV reproduces the observed data remarkably well. In 
this dispersion relation, the averaged excitation energy, the band width of each excitation mode, and the energy 
difference between the acoustic and optic modes are related with J, J cp -and J a , respectively. J and J c obtained in 
this experiment are consistent with those obtained by Eccleston-et-aLt£l The sign of J a is consistent with that in 
Lai4_ x Ca x Cu2404x which shows a long range magnetic ordering £ir£!3 
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RPA treatment gives the dynamic structure factor0 as the following 



S(Q, oj) oc F 2 (Q) ■ R(T) ■ J ■ [1 - cos(Q • R) 



[1 + COs(p • T + (j))]- 



.acoustic 



(q) 



2 (<z)) 



[i - cos( P • T + W^^s (« - - optic (g)) 



(3) 



where cj> = 2wl8, Q=q + t (t is a reciprocal wave vector), and i? is the vector connecting individual spins within 
a dimer. The factor [1 — cos(Q • R)] can be seen in the structure factor for isolated dimers. The acoustic and optic 
modes can be distinguished by the factor [1 ± cos(p • r + <fi)]. The solid and broken lines in Figs. 3(b)-6(b) represent 
the energy-integrated intensities of the acoustic and optic modes calculated using Eq. (3), respectively. J, J a , and X, 
are fixed at the values mentioned above. 5 is fixed at 0.77 which was determined by x-ray measurements at 50 Kx3 
Only adjustable parameter is the scale factor. The same scale factor was used for the intensities at (2,0,L), (3,0, L), 
(H,0, —0.2), and (H,0, —0.7). The energy-integrated intensities are described with the RPA theory reasonably well. 



3. Temperature dependence 

Figure 8 shows constant-Q scans at (2,0,-0.1) and (2,0,-0.6) as a function of temperature. The solid lines 
are the fits to two Gaussians. In the fitting the higher excitation energy at (2,0,-0.1) and the lower excitation 
energy at (2, 0, —0.6) are fixed, which is confirmed with more intense excitation peaks at (3, 0, —0.1) and (3, 0, —0.6). 
It was assumed that the peak width and the intensity ratio of the two excitations are temperature independent. 
Temperature dependence of excitation energies and energy-integrated intensities are plotted in Fig. 9. The lower 
excitation energy at (2,0,-0.1) becomes higher and the higher excitation energy at (2,0,-0.6) becomes lower with 
increasing temperature. Furthermore, the intensities decrease with increasing temperature. Since there is a thermal 
factor R(T) in Eqs. (2) and (3), the RPA calculation predicts the temperature dependence of the dispersion relation 
and the energy- integrated intensities. The dispersion becomes flatter with increasing temperature as shown in the 
inset of Fig. 9. because the interdimer couplings become negligible due to large thermal fluctuations and the dimers 
behave like isolated dimers. The intensity is proportional to R(T) which decreases with increasing temperature. The 
solid lines in Fig. 9 are the results of calculations using Eq. (2) and (3). All the parameters J, J a , J c , 5, and 
scale factor are fixed at the values as determined above. The temperature dependences are fitted with the RPA 
calculation reasonably well although the observed intensity decreases more quickly than the calculated intensity does. 
This indicates that other factors in addition to thermal fluctuations disturb the dimers and the ordering of Cu 2+ 
and ZR singlet- A gradual destruction of the ordering with increasing temperature suggested by NMRH and x-ray 
measurementsEj might cause the behavior. 



4- Discussion 



Now we_discuss the relation between magnetic dimerization and structural distortion. As mentioned in Sec. I, 
Cox et ciIIlB performed synchrotron x-ray measurements on Sri4Cu2404i and observed superlattice Bragg peaks at 
(0,0, Z/4) and (0,0, 1/2) which|-o*igmate from a structural distortion in the chain. This lattice distortion is probably 
related with the NMR resultsaO since the superlattice intensity also decreases gradually with increasing temperature. 
It should be pointed out that the x-ray experiments properly detected the coupling along the a axis but there is some 
disagreement along the c axis. As mentioned above, the structural distortion has correlations with two and four times 
periodicities along the c axis although magnetic dimers have correlations with five times periodicity^- ,A discrepancy 
between magnetic and structural correlations was also reported in the spin-Peierls state in CuGe03.c2l The magnetic 
excitations have a minimum energy at (0,1,^) but a superlattice reflection at (i,l,|). An important point is that the 
x-ray results are for the ground state and neutron results are for the excited state. 

We have shown that the dimer model shown in Fig. 7 explains the magnetic excitations reasonably well. The 
dimerized state is ascribed to the ordering of Cu 2+ and ZR singlet which originates from an ordering of the localized 
holes. Hole ordering or charge ordering occurs in various 3d transition metal oxides. A quasi-two-dimensional hole 
ordering of the doped holes was reported in L^NiCU+tS" and La2-a-3/Nd x Srj / CuC , 4EM. In the phase a static ordering 
of antiferromagnetic stripe is separated by charge ordered domain walls. In the latter compound the ordered phase 
can be described as a stripe ordering of CuS" and ZR singlet. Manganese oxides, fop-example Lai-^Ca^MnOs, also 
show a charge ordering of Mn 3+ and Mn 4+ H3 In this system spin, charge, and orbitalEU degrees of freedom are closely 
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related, which results in interesting phenomena such as colossal magnetoresistance.Ea NaV20s, which consists oL two- 
leg ladders of vanadium ions, has been studied extensively since it shows a sin skit ground state below 34 K.El The 
phase transition was first considered to be a conventional spin-Peierls transition.^ But the detailed inelastic neutron 
scattering experiment using single crystal sample revealed that the. magnetic excitations in the dimerizcd phase could 
not be understood by a simple one-jdimensional dimerized model.! 41 ! NMR experiment also indicated a charge ordering 
of V 4+ and V 5+ at the transitions Stimulated by these experimental findings, new approaches to understand the 
origin of the phase transition have been made theoreticallyOc-l It is claimed that the transition is indeed closely 
related with a charge ordering of V 4+ and V 5+ ions. These various properties suggest that an intimate connection 
between hole/charge ordering and magnetic ordering is a common feature in the strongly correlated 3c? transition 
metal oxides. Theoretical studies to explain the hole ordering and the dimerized state in Sri4Cu2404i are highly 
desirable. 



B. Substitution effect in Sri4_ :c Ca a ,Cu2404i (x=3) 

We have also studied the substitution effect on the dimerized state. The dimerized state in the chain is most stable 
in pure Sri4Cu2404i. When Y is substituted for Sr, the number of the holes is decreased and the number of the Cu 2+ 
spins is increased. It was reported that Y substitution makes the magnetic excitation peaks broaderXJ As mentioned 
in Sec. I, the Ca substitution for Sr also reduces the number of the holes in the chain gradually-Enl 1 ] Therefore, the 
Ca substitution is expected to affect the dimerized state in the chain. 

Figure 10 shows typical inelastic neutron spectra at (3,0,L) measured at 8 K in SrnCa3Cu240«, The magnetic 
excitation peaks become broader with Ca substitution, which is similar to the case of Y substitution.LJ It is noted that 
two excitation branches are hardly resolved. However, the excitation energies are similar to those in Sri4Cu24C>4i. 
The solid lines at (3, 0, —0.15), (3, 0, —0.30), and (3, 0, —0.35) are fits to two Gaussians and that at (3, 0, —0.75) is a 
fit to a single Gaussian by assuming that the dispersion relation is the same as that in Sri4Cu2404i. The spectra are 
reasonably described with the simple model. As mentioned in Sec. III-A, the averaged excitation energy, the band 
width of each excitation mode, and the energy difference between the acoustic and optic modes are related with J, 
J c , and J a , respectively. This indicates that the coupling constants J, J c , and J a are almost unchanged with Ca 
substitution although the dimerized state becomes unstable. This behavior would be explained as the following. The 
hole ordering sensitively depends on the hole number. Accordingly, the long-range dimer formation becomes disturbed 
with Ca or Y substitution. nr-\ 

In Sri4_ a Y a; Cu2404i system, magnetic excitations are remarkably affected with Y substitution. EjO In the x=0.25 
sample the excitation peaks become broader although the dispersion relation is almost unchanged. This result suggests 
that J, J c , and J a are almost unchanged although the hole ordering is disturbed with a small amount of Y. In the 
x=l sample the excitation peaks become much broader, suggesting that the hole ordering becomes more unstable. 
An interesting behavior is that the averaged energy of the two excitations becomes reduced to ~9 meV, indicating 
that J is decreased in Sri3YiCu2404i. The changes of the bandwidth of each excited state, the difference in energy 
between two excited state, and the periodicity of the dispersion relation are difficult to be determined because the 
peaks are too broad to be resolved as two excitations. A puzzling feature is that J is decreased in Sri 3 YiCue4.04i 
although the lattice constant c, which affects the exchange constant, is almost independent of Y concentration!^ 

When the hole s are re moved further, the hole ordering cannot be observed anymore and a long-range magnetic 
ordering appears. EirEla The Cu02 chains in the hole removed Sri4Cu2404i can be considered as ferromagnetic spin 
chains which are weakly coupled with antiferromagnetic interchain coupling along the a and b axes. 

In conclusion, neutron scattering experiments have been performed on Sri4Cu2404i in order to study the dimerized 
state in the CuC-2 chains. The uj-Q dispersion relation as well as the integrated intensities can be reasonably described 
by RPA calculation with J = ll meV, </ c =0.75 meV, and J a =0.75 meV. The dimer configuration indicates a quasi-two- 
dimensional hole ordering, which results in an ordering of magnetic Cu 2+ with S = \ and nonmagnetic ZR singlet. 
The hole ordering and the dimerized state are most stable in Sri4Cu2404i. The Ca substitution for Sr sites makes 
the excitation peak broader probably because the hole ordering becomes unstable and the long-range dimer formation 
becomes disturbed. 
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FIG. 1. Structure of the CuCb chains (a) and the CU2O3 ladders (b) in Sri4Cu2404i. The lower part of the figure shows 
two dimer models as described in Section I: (c) Model I and (d) Model II. The squares denote Zhang-Rice singlet sites. Figure 
from ref. 18. 

FIG. 2. Typical inelastic neutron spectra at (H ,0,L) at 15 K in Sri4Cu24 04i. The solid lines are the results of fits to a single 
Gaussian or two Gaussians. 

FIG. 3. Observed and calculated energies (a) and intensities (b) at (2,0,L) measured at 15 K in Sri4Cu2404i. The solid and 
broken lines, which are fits to Eqs. (2) and (3), represent the acoustic and optic modes, respectively. 

FIG. 4. Observed and calculated energies (a) and intensities (b) at (3,0,L) measured at 15 K in Sri4Cu2404i. The solid and 
broken lines, which are fits to Eqs. (2) and (3), represent the acoustic and optic modes, respectively. 

FIG. 5. Observed and calculated energies (a) and intensities (b) at (H, 0, —0.2) measured at 15 K in Sri4Cu2404i. The solid 
and broken lines, which are fits to Eqs. (2) and (3), represent the acoustic and optic modes, respectively. 

FIG. 6. Observed and calculated energies (a) and intensities (b) at (H, 0, —0.7) measured at 15 K in Sri4Cu2404i. The solid 
and broken lines, which are fits to Eqs. (2) and (3), represent the acoustic and optic modes, respectively. 

FIG. 7. A proposed model for the dimerized state and the ordering of Cu 2+ and ZR singlet in the ac plane. 

FIG. 8. Inelastic neutron spectra at (2,0,-0.1) (T=15 and 75 K) and (2,0,-0.6) (T=15 and 100 K) in Sri 4 Cu240 4 i. The 
solid and broken lines are the results of fits to two Gaussians. 

FIG. 9. Temperature dependence of the energy-integrated intensity and energy of the gap excitations at (2,0,-0.1) and 
(2,0, —0.6) in Sri4Cu24 04i. The solid lines are fits to Eqs. (2) and (3). The inset shows the calculated dispersion at 15 and 
100 K. 

FIG. 10. Typical inelastic neutron spectra at (3,0,-L) at 8 K in SrnCa3Cu2404i. The solid lines are the results of fits to a 
single Gaussian or two Gaussians. 
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